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We report the first lattice QCD calculation using the almost physical pion mass — 149 MeV 
that agrees with experiment for four fundamental isovector observables characterizing the gross 
structure of the nucleon: the Dirac and Pauli radii, the magnetic moment, and the quark momentum 
fraction. The key to this success is excluding the contributions of excited states. An analogous 
calculation of the nucleon axial charge governing beta decay fails to agree with experiment, and we 
discuss possible sources of error. 

PACS numbers: 12.38. Gc, 13.60.Fz 



INTRODUCTION 

Lattice QCD is the only known rigorous framework 
for ab-initio calculation of the structure of protons and 
neutrons with controllable errors. It can provide quanti- 
tative answers to both fundamental questions such as the 
quark and gluon composition of the nucleon spin and phe- 
nomenological questions such as the sensitivity of modern 
detectors to physics beyond the Standard Model (BSM), 
to fundamental symmetry violations, and to hypothetical 
dark matter particles [TH3]- However, with current com- 
puter resources, its predictive power is limited by uncer- 
tainties arising from heavier than physical quark masses, 
finite lattice spacing and volume, incomplete removal of 
excited states, and omission of disconnected contractions. 
Therefore, until exhaustive lattice calculations remove 
these uncertainties, reproducing several well-known ex- 
perimental observables is an important way to increase 
confidence in lattice QCD predictions. 

Significant effort has been focused on lattice calcula- 
tions of several isovector quantified] such as the Dirac 
and Pauli radii (r 2 2 Y t the axial charge gA, and the quark 
momentum fraction (x) u —d\ 

(j/\rfq\p) = V [FliQ 2 )^ + F«(Q 2 ) l -^]u p , 

2M (1) 

FUQ 2 ) Q2 ~ *?,a(o)(i - - 6 {rh) q Q 2 + o(Q 4 )), 

(pWli^D^p) = (x) q Up-y^p^Up, (2) 

(p\qi"i 5 q\p) = 9a u P Yi 5 u p , (3) 
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1 So-called disconnected contractions, which are expensive to com- 
pute, cancel in isovector observables, making them ideal for ver- 
ifying lattice QCD. 



where Q 2 = —q 2 = —{p' — p) 2 and u p , u p r are nucleon 
spinors. Although some success has been achieved [4]- 
E] , past results rely heavily on large extrapolations using 
Chiral Perturbation Theory (ChPT) yielding potentially 
uncontrollable corrections. This is particularly problem- 
atic for some observables, e.g., (rf t2 ) v and (x) u -d, for 
which ChPT predicts rapid change towards the chiral 
regime, making extrapolations very difficult. For ex- 
ample, in typical lattice calculations with pion masses 
> 250 MeV, prior to extrapolation to m* 8 « 135 MeV, 
(r\) v is underestimated by « 50% [6"H5I[TU]. (x) u -d over- 
estimated by 30 — 60% [3 HU UH , and qa underestimated 
by sa 10% [HQS [14] , compared with experiment. These 
glaring discrepancies and the dependence on large extrap- 
olations clearly indicate the need for calculations near the 
physical pion mass. Moreover, it has recently been found 
that excited-state effects become worse with decreasing 
pion mass [15], and their careful analysis is required be- 
fore even attempting extrapolations in the pion mass to- 
wards the physical point using ChPT. 



In this paper, we report the first lattice QCD calcu- 
lation of nucleon structure using pion masses as light as 
TOtt = 149 MeV and thus very close to the physical value; 
therefore, our results rely much less on ChPT extrapo- 
lations than previous calculations. For each ensemble, 
we remove excited-state contaminations by varying the 
source-sink separation in the range « 0.9 ... 1.4 fm and 
apply the summation method |16] to extract the ground 
state matrix elements. We observe remarkable agreement 
with experiment of the isovector Dirac and Pauli radii, 
the anomalous magnetic moment, and the quark momen- 
tum fraction, all computed with the same methodology. 
However, as we will see later, the axial charge gA is still 
underpredicted and requires further studies. 
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LATTICE RESULTS 



We perform calculations using ten ensembles of lat- 
tice QCD gauge fields generated with C(a 2 )-improved 
Symanzik gauge action and tree-level clover-improved 
Wilson fermion action using 2-HEX stout gauge 
links |17j . Two light u and d (with m u = m^) and one 
heavier strange (m a ^> m u> d, tuned to be close to physi- 
cal) quark flavors are simulated fully dynamically, while 
effects of heavier quarks are neglected. In addition to 
varying the pion mass in the range 149 . . . 357 MeV, we 
include different spatial volumes, time extents of the lat- 
tice, and one ensemble with a smaller lattice spacing in 
order to estimate the size of the corresponding systematic 
effects; see Tab.[TJ Nucleon matrix elements are extracted 
from nucleon 3-point functions that are computed with 
the standard sequential source method. Nucleon field 
operators are optimized to overlap as much as possible 
with the single-nucleon ground state at rest by tuning the 
spatial width of Gaussian smeared quark sources. This 
width is kept approximately constant at different pion 
masses to avoid possible bias due to different smearing. 
In order to discriminate between the ground and excited- 
state matrix elements on a Euclidean lattice, we vary the 
timelike distance At between nucleon sources and sinks in 
the 3-point functions. With increasing At, excited states 
in 3-point functions are suppressed as ~ e -A_EAt/2 anc j 
disappear in the At — > oo limit, where AE is the (po- 
tentially m^-dependent) energy gap to the closest con- 
tributing state. However, using a large source-sink sepa- 
ration is impractical, since statistical noise grows rapidly 
with At. Instead, we combine calculations with three 
values of At w 0.9, 1.2, 1.4 fm using the summation 
method |16j . which benefits from improved asymptotic 
behavior [TSl [H] and which we find the most reliable and 
robust method with the presently existing statistics. In a 
separate publication [20] we report detailed studies and 
comparisons of different methods, including GPoF [21] 
and n-state fits. 

We present and discuss our results for all the observ- 
ables in Figs. [Tfj5] below. To emphasize the importance 
of controlling excited states, in the same figures we also 
show the standard plateau method results with fixed At 
(open symbols) for the three lightest used in fits. As 
this separation is increased, the data points consistently 
approach our final results, while their error bars increase 
as expected. As illustrated by the data, the effect of the 
excited states may be very dramatic, especially for (rf) v , 
K v {rl) v and (x) u - d - 

Our lightest pion mass value — 149(1) MeV is only 
w 10% higher than the physical pion mass in the isospin 
limit m-n- ps 134.8 MeV [22]. Since we still have to extrap- 
olate our data in m w , we select two additional ensembles 
with the next-lightest pion masses 202(1), 254(1) MeV 
and the largest spatial volumes to keep model and sys- 
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m 3 is tuned to be close to physical. 
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TABLE II. ChPT extrapolations to the isospin-limit m^ hys = 
134.8 MeV [22] using data points 149 < < 254 MeV. The 
footnotes describe details of ChPT and fixed parameters in 
the fits. The right column shows the deviation from experi- 
ment, relative to the lattice uncertainty. 
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a SSE C(e 3 ) [26] with fixed F° = 86.2 MeV [27], 
A = 293 MeV [ID], g° A = 1.26 lil l28ti3P1 , c A = 1.5 HDI . 

b No A resonance (c A = 0). 

c Includes (NA) M l transition with c v = -2.5 GeV -1 [31] 

d Includes higher-order "core" term I32| . 

e BChPT 0(p 2 ) [33] with fixed F° = 86.2 MeV [27J, 
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f SSE 0(e 3 ) [28] with fixed F° = 86.2 MeV g7J, 
A = 293 MeV QDJ, c A = 1.5 QDJ, 91 = 2.5 [35]. 
s From r p E [23] and (r|) n [24]. Using Ref. [24] for both (r|)P> n 
results in (rf) v = 0.640(9) (higher exp. point in Fig.JTJ. 



tematic bias minimal. Table [IT] contains details of fits, 
results and comparison with experiment. To check the 
consistency of the low-energy theory, we repeat chiral fits 
using the full available range of = 149 . . . 357 MeV. 
In general, the results from < 350 MeV and < 
250 MeV fits agree well within error bars, demonstrating 
good convergence of ChPT. 

To explore nucleon electromagnetic structure, we com- 
pute matrix elements of the quark-nonsinglet vector cur- 
rent (p'\ [uj^u — dj^d] \p) between polarized proton (uud) 
states with different momenta p, p'. We extract isovector 
Dirac and Pauli form factors F^ 2 {Q 2 ) from these matrix 



elements, and then the "radii" (r\ 
using dipole fits F 12 (Q 2 ) 
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FIG. 1. Isovector Dirac radius (r\) v . Fits to the blue solid 
square and diamond points are described in Tab. [IT] and the 
same fits applied to the full set of solid points are shown for 
comparison. The two experimental points are from PDG |24| 
and the fip Lamb shift [23] . The series of open symbols show 
data before the removal of excited states, with fixed source- 
sink separation At increasing from right to left. Their error 
bars reflect only statistical errors, which grow with At. 
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FIG. 2. Isovector anomalous magnetic moment k v . See 
caption of Fig. [JJ 



< Q 2 < 0.5 GeV 2 . 

The Dirac radius {r\) v is shown in Fig.JT] We compare 
it to the experimental value (r 2 ) v = (rf) p — (r\) n , where 
(r 2 ) p ' n = (r%) p ' n - , with the error bar dominated 

by the uncertainty in (r 2 E ) p , the proton electric charge 
radius. We show two experimental values for (r 2 )" in 
Fig. [JJ which correspond to two inconsistent values for 
(r%) p : the PDG value [24] and the recent and controver- 
sial result from measurement of the up Lamb shift [33J . 
Relative to the lattice uncertainty, the extrapolated value 
deviates from the /j,p Lamb shift value by —0.07(7 and 



FIG. 3. Isovector Pauli radius k v ' (r^Y ■ See caption of Fig. [I] 
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FIG. 4. Isovector quark momentum fraction {x) u -d- See 
caption of Fig. [JJ 



from the PDG value by — 2.17a. In addition, we check 
the low-energy QCD dynamics by repeating the ChPT 
fit without the A-resonance, shown as the dotted lines 
in Fig. [JJ In this case, we observe somewhat worse fit 
quality (see Tab. |H{ line 2), especially when the full 
range is included in the fit, demonstrating the relevance 
of the A-resonance. 

In a similar fashion, we extract the isovector anoma- 
lous magnetic moment K v — k p — n n and the Pauli radius 
(r 2 ) v = (K p (r%) p -K n (r%) n )/(K p -K n ) from the Pauli form 
factor F2(Q 2 ). In this case the results are less precise 
because the forward values F 2 (0) and are extrapo- 



lated using the dipole form F£ (Q ) 



(i+QyM% 2 ) 



,-. Since 



2 ^ J_ 
min L 2 ' 



the minimal value Q 2 > scales roughly as Q 
the Q 2 fits are less precise on lattices with smaller spa- 
tial volumes. This explains the significant increase of 
error bars in Figs. [2J [3] going from 32 3 to 24 3 lattices 
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FIG. 5. Axial charge g A - See caption of Fig. [T] 

(e.g., at 77i fa 250 MeV), compared to the correspond- 
ing error bars of {r\) v in Fig. [l] However, our results 
at the lightest pion mass are computed with the largest 
V 3 ps (5.6 fm) 3 yielding the smallest Q^ in ps 0.05 GeV 2 
and are the least prone to this problem. For simplic- 
ity, we fit and compare to experiment the combination 
/^(rf)", because this quantity is more natural for ChPT 
than the radius (r|)" itself. We achieve remarkable agree- 
ment both for the anomalous magnetic moment n v and 
the isovector Pauli radius (r%) v . 

We compute the isovector quark momentum fraction^] 
{x)u-d from the forward matrix element of the operator 
in Eq.[2]and renormalize the lattice value to the standard 
MS(2 GeV) scheme using the RI'/MOM method [36]. 
In Fig. [4] we show ChPT extrapolation curves together 
with the CTEQ6 phenomenological value [25 . Achieving 
agreement between the lattice result and phenomenology 
for (x) u —d is one of the most important accomplishments 
of this paper: previous lattice calculations [12l ESJ [37] 
consistently overestimated the phenomenological value. 
Our removal of excited state contamination eliminated 
the discrepancy. 

In stark contrast with the observables discussed above, 
the computed value of the nucleon axial charge g A 
strongly disagrees with g c * p = 1.2701(25) [53]. There 
is a notable tendency that g 1 ^ monotonically decreases 
with TOtt — > m phys for < 250 MeV, increasing the 
discrepancy with g A * p ■ It is remarkable that our re- 
sults for 777^ > 250 MeV agree with the recent calcu- 
lation [5] that used similar techniques, but decrease dra- 
matically for lighter pion masses. Apparently, there is 



2 (a:) u — d = {x) u — is understood as the momentum fraction 
carried by quarks and antiquarks, i.e. (x) q = dx x (f q (x) + 
fq(x)), where fqt fi is a parton distribution function (PDF). 



a source of bias at light pion masses that affects the 
axial charge significantly more than the other observ- 
ables discussed above. According to Ref. [13], finite 
volume effects (FVE) may lead to a 9% bias in g A at 
m^Ls ps 4.5 in calculations with domain wall fcrmions, 
and as much as ps 25% in calculations with Nf = 2 
flavors of Wilson fermions. However, from our data, 
it seems very unlikely that FVE can explain the de- 
viation of our near-physical pion mass result from ex- 
periment: assuming that the FVE correction scales as 
S FYE g A ~ e~ m - L analogously to Ref. [JJ] and using the 
gA values from the two lattices with ps 250 MeV that 
differ only in spatial volume [m^Ls = 3.6 and 4.8), we 
obtain an estimate 5 FVE g A = 0.02(8) for m n = 149 MeV, 
L s = 5.6 fm, which is much smaller than the observed 
discrepancy Sg A ~ —0.3. On the other hand, the data 
with 771^ < 250 MeV hint that thermal states [35] may 
cause a significant bias. For example, pion states are sup- 
pressed only by a factor e r m ^/ T = e -m v L t ^ anc j ma y ^ e 
natural candidates to make large contributions to matrix 
elements of the axial current. Although the statistics are 
not sufficient to draw a solid conclusion, at given pion 
mass the central values of g A in Fig. [5] monotonically rise 
with m^Lt = 3.6. ..14.3 (see Tab. ffl, supporting this hy- 
pothesis. This effect appears to weaken in the vicinity of 
m^L t ~ 7.2. 

DISCUSSION 

The crucial departure of this work from previous calcu- 
lations lies in the careful control of excited-state contam- 
inations in nucleon matrix elements on lattices almost 
all the way down to the physical pion mass. We have 
shown that, although this is multiplicatively more com- 
putationally expensive (in this case, x3), such control is 
essential for extracting correct results from lattice QCD 
at the physical pion mass and reproducing experimental 
values of nucleon structure observables. 

We have carried out calculations with a range of pion 
masses and consistently applied the robust summation 
technique to remove excited states. In addition, we have 
performed some limited checks of finite volume and dis- 
cretization effects. For the first time, we have achieved 
agreement with experiment for all the reported quanti- 
ties except the axial charge g A . This agreement is strong 
evidence that this combination of methods is sufficient to 
correctly estimate and remove systematic errors in most 
nucleon structure calculations. Our results suggest that 
the disagreement in the axial charge value is due to the 
insufficient temporal extent of the lattice, which is equiv- 
alent to having non-zero temperature in the system. 

Additional statistics are necessary to improve the pre- 
cision of the reported results, since the errors are still 
much larger than in experiments. Nevertheless, the re- 
markable precision of the computed isovector Dirac ra- 
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dius {r\) v is already comparable to the discrepancy be- 
tween the two contradictory experimental values, and fu- 
ture improvements may be able to contribute to a reso- 
lution of this discrepancy. 

More lattice spacings and further volumes and tempo- 
ral extents at pion masses below 200 MeV are required 
to put stronger bounds on other sources of systematic 
uncertainty. With appropriately increased statistics, our 
control over excited-state contamination could be cross- 
checked by further varying the nucleon source-sink sepa- 
ration and by comparing against other excited state re- 
moval techniques. 

It is encouraging that Lattice QCD, for many impor- 
tant benchmark quantities, now is in good agreement 
with experimental results. 
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